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A  TRAFFIC  COUNTING  DISTRIBUTION 


0 .  Introduction 

Several  authors^*5^  have  suggested  the  use  of  traffic  counting 
distributions  in  which  a  minimum  spacing  or  headway  between  units  of 
traffic  is  taken  into  account.  Observers  of  medium  and  high-density 
traffic  flows  along  highways,  freeways  and  air  lanes  have  evidence  to 
support  such  a  claim.  In  the  case  of  road  traffic  a  minimum  spacing 
may  be  governed  by  the  finite  car  sizes,  the  minimum  headways  by  the 
velocity  and  spacing  characteristics  of  drivers.  For  reasons  of  safety 
airplanes  may  have  to  be  separated  by  a  minimum  space  or  time  interval. 

While  at  least  one  author^  has  pointed  out  that  the  relative  merits 
of  minimum  spacings  are  evident  in  the  measurement  of  inter -vehicle 
distances,  this  may  not  be  the  case  in  measurements  of  headways,  i.  e.  , 
the  time  between  the  passing  of  two  vehicles.  Leaving  aside,  momen¬ 
tarily,  the  question  of  the  domain  of  definition  of  the  velocity  distribution 
and  its  fundamental  effect  on  space-time  measurements,  the  author  of 
this  paper  feels  that  the  probability  distribution  of  inter- vehicle  headways 
which  includes  a  minimum  headway  may  serve  as  a  valuable  approximation 
to  those  cases  where  small  headways  are  merely  improbable. 

A  number  of  probability  distributions  of  inter-vehicle  spacings  or 
headways  have  been  proposed;  in  many  cases  the  corresponding  counting 
distributions  have  been  found.  When  the  inter-vehicle  spacing  is  expo¬ 
nentially  distributed,  the  discrete  counting  distribution  is  the  well  known 
Poisson  case.  When  the  probability  density  distribution  of  inter- vehicle 
spacings  is  Erlang,  Gamma  or  Type  III,  the  corresponding  counting 
distribution  is  that  of  the  Generalized  Poisson  (one  starts  to  count  with 
the  passing  of  a  vehicle)  or  the  state  probabilities  discussed  in  a  ps.per 
by  Jewell^*  (one  starts  to  count  at  a  random  instant  of  time),  '.n  this 
paper  the  author  obtains  solutions  for  the  discrete  probability  distribu¬ 
tions  and  the  average  and  variance  of  vehicle  counts  in  an  interval  when 
the  inter- vehicle  spacings  have  the  probability  density  distribution  of 
Equation  ( 1) . 


The  State  Probabilitie s 

Let  the  probability  distribution  of  the  spacing  between  units  of 
traffic  be  given  by 


a(x)  =  0 

0  <  x  <  Xq 

(la) 

=  He^-V 

*0  1  x 

(lb) 

i.e.  ,  the  probability  than  an  event  occurs  in  dx  is  pdx  provided 

x  is 

greater  than  a  minimum  spacing  Xq. 

By  standard  arguments^) 

which 

neglect  terms  proportional  to  square  or  higher  powers  of  dx  we 
that  the  Chapman- Kolmogorov  equations  for  the  state  probability 

find 

Pn(x) 

of  a  count  of  n  events  in  an  interval 

,  * 

x  is  given  by 

dP0 

~ST  =  1 

o 

1  A 

X 

1  A 

X 

O 

(2a) 

=  *PP0 

xo  5  x 

(2b) 

dp 

n  _ 

~ar  -  0 

0  <  x  <  nxQ 

(3a) 

=  MVi-l 

nx0  -  X1  (n  +  1)xo 

1 3b) 

=  MVl-M*  n 

(n  +  1)xq  <  x 

(3c) 

In  solving  these  equations  of  the  birth  process  one  can  obtain  the 
constants  of  integration  by  making  use  of  the  condition  that  the  sum  of 
the  state  probabilities  add  to  one.  For  a  given  interval  x  and  minimum 
spacing,  Xq  ,  the  sum  extends  over  a  finite  number  (say,  n  =  N)  which 

Pn(x)  corresponds  to  An<x)  in  Reference  6;  i.e.,  the  counting  origin 
starts  with  the  passing  of  a  vehicle. 


-2- 


represents  the  fact  that  a  maximum  number  of  events  can  occur  in  the 
interval  x.  N  is  just  the  integral  part  of  xXq  *  which  we  denote  by  square 
brackets, 


N 


xx. 


-1 


Hence, 


N 

n=0 


(4a) 


(4b) 


Although  the  algebra  becomes  involved  for  n  >  1,  the  solutions  of  Equations 
(2),  (3)  are  obtained  in  a  straight-forward  manner .  The  solution  of  Pq(x)  i® 

1  in  the  interval  (O.Xq)  and  e"M-(x_xo)  for  x  >  xQ.  This  exponential  term  is 
just  the  probability  that  the  spacing  between  units  of  traffic  is  greater  than  x, 
i.  e.  ,  the  tail  of  the  distribution  of  a(x)  in  Equation  (1). 

The  probability  of  one  event  in  (O.Xq)  is  zero  since  units  of  traffic  must 
be  separated  by  at  least  Xq.  The  probability  of  one  unit  in  x(x^<  x  <  2xQ) 
is  just  1  -  e  X°^  since  there  can  only  be  zero  or  one  unit  in  an  interval 

less  than  2x^  in  length.  The  probability  of  one  unit  in  x(2xq  <  x)  is  the 
solution  of  (3b)  when  one  substitutes  1  for  n  and  e  *^X  X°^  for  Pq(x).  In 
summary  we  find  that 

p^y)  0  5  y  5  xo  *5a-' 

Pi(y)  =  1  -  e"P’(y_xo^  x0f  y  <  2xq  (5b) 

=  p(y  -  2xq)  e’^'2^  +  e-P(y2xo) 

-e'^y  ■  Xq)  2x0  <  y  (5c) 

In  general  we  find  that  the  probability  distribution  is  composed  of  three 
parts- 
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(6a) 


pn(y)  =  0 


0  <  y  <  nx. 


=  y  f  n;  p(y  -  nxQ)  /(n-1)!  nxQ  <y<(n+l)xc 

=  y  n;  p(y -nxQ)l  /  (n  -  1)!  (n+l)xQ<y 


-V  n  +  1; 


p(y-(n+l)x 


0)  ]  /  ni 


where  y(n;  x)  is  the  incomplete  Gamma  function  of  order  n  and  argument  x' 
The  discrete  distribution  Pn(y)  is  continuous  over  the  region  0  <  y  <  oo  and 
is  plotted  in  Figure  1  for  several  values  of  N  =  [yx^  *]  and  average  count 


W.  Feller'  has  pointed  out  that  the  probability  distribution  Pn(y)  can 

be  expressed  as  the  difference  between  two  expressions  F^  ^(y)  and  F^(y) 

which  are,  respectively,  the  probability  that  the  interval  between  the  n 
s  t 

and  (n+1)  event  is  less  than  or  equal  to  y.  Since  we  are  considering 
counting  intervals  which  begin  with  the  passing  of  a  traffic  unit,  Feller's 
expression  for  the  probability  that  the  time  to  the  next  event  is  less  than  or 
equal  to  x  becomes 


Fq(x)  =  0 


=  1  -  e-»1<x"xo) 


0  5  x  1  xo 


x  >  xo 


and  the  probability  that  the  time  between  the  passing  of  the  n**1  and  (n  +  l)8* 
unit  is  less  than  or  equal  to  y  is  just 


In  this  expression  and  the  ones  which  immediately  follow,  Feller  uses 
F^x)  =  1  -  e  x  >  0,  i.e.  ,  a  counting  origin  randomly  placed  outside 
the  minimum  separation  interval  x«. 
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Counts,  n 

THE  PROBABILITY  OF  n  COUNTS  IN  s  AS  A  FUNCTION  OF  n 


ry 

Fn<y>  =  )  Fn-l<y-x)dF0(x) 

‘o 

=  ^  KFn  l(y  -  x)e‘^X 


dx 


(8a) 


(8b) 


(8c) 


=  y  [n  +  1;  p(y-(n+l)xQ)]  /  n! 

_  i 

The  difference,  Fr  ^(y)  -  FR(y)  is  just  Equation  (ob). 


Average  and  Variance  of  Traffic  Counts 

Standard  renewal  arguments  can  be  invoked  to  derive  the  average 
and  variance  of  the  distribution  of  Equation  (6).  They  can  also  be 
derived  by  inversion  of  the  derivatives  of  the  double  transform 


00  00 


P(s;z)  =  ^  J  e"8y  z11  Pn(y)  dy 


(9) 


n=0  0 

Fortunately/"^  ^  the  expression  for  P(s;  z)  and  moments  can  be 
obtained  in  terms  of  the  Laplace  transform  a(s)  of  the  inter-event 
distribution,  a(x),  in  Equation  (1): 


a(s)  =  ^  e  8X  a(x)dx  = 


(10) 


The  expression  for  P(s;  z)  is  then 


** 


In  the  notation  of  incomplete  Poisson  sums,  -^^-’x^ |  =  j(x)> 

j=n 

the  expressions  of  Equations  (6)  and  (8c)  are  of  the  form  1  -  En(py-(n+l)xQ). 

On  expanding  P(s;z)  in  powers  of  z  we  obtain  the  Laplace  transform  of  p  „(y). 
with  inver  se  transform  again  equal  to  the  solutions  of  Equation  (6). 
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J. 

s 


(ID 


(  e+y.  -pc~BX °  \ 

^  s  +  jj.  -  pze  o  ' 

and  the  Laplace  transform  of  the  mean  count,  n(s)  is  just 

^(s)  =  a(s)  [  s  -  s”a(s)]_1 

•SXnr  2  ,  -  6X„i  -  1 

=  pe  °[s  +  (jls  -  p.se  °J 

with  inverse  transform  equal  to  the  average  count  in  y 

V  Vt JSM-Cy  -  jx  )] 
n(y)  -  2  - 0^-m - 

j  =  l 

By  similar  arguments  the  variance  is  found  to  be 
N 

Var  (n(y) )  =  ^  (fjljl  V  [  j iH-(y- jxQ)] 

j-1 


(12a) 

(12b) 


(13a) 


(13b) 


Figure  2  is  a  plot  of  the  variance  versus  the  mean  for  values  of 
N  =  5,  10,  20.  The  maximum  average  counts  are  4,  9,  19.  This  maximum, 
N-l,  is  due  to  our  choice  of  origins  where  one  unit  of  traffic  passes  just 
before  time  zero.  The  distribution,  pn(x),  describes  the  probability  of  n 
events  in  addition  to  the  one  counted  at  time  zero;  hence,  the  maximum 

value  ot  rv(x)  reflects  this  feature  of  being  one  less  than  the  integral  part 

,  -1 

of  XX 

o 

Several  interesting  properties  of  the  variance  to  mean  ratio  can  be 
observed,  ui  t.^e  first  place  we  note  that  for  low  average  counts  in  an 
interval  x  the  limiting  variance  to  mean  ratio  is  unity,  i.e.  ,  that  of  the 
Poisson  distribution.  When  average  spacings  between  cars  are  much  larger 
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Fig.2-THE  VARIANCE  VERSUS  THE  MEAN. 


than  the  minimum  spacing  xq  the  low-density  Poisson  counts  are 
obtained.  As  the  mean  count  increases  the  variance  also  increases, 
reaches  a  maximum  and  then  decreases  to  zero.  Var  n(x)  is  always 
less  than  n(x)  and  the  probability  distribution  pn(x)  has  a  "maximum¬ 
packing"  property  in  the  sense  that  the  variance  of  traffic  counts  is 
zero  when  the  mean  count  in  x  equals  the  maximum  count  N  =  [  xxq"  1  ]  . 
As  we  have  already  seen  in  Equation  (6)  pn(x)  =  0  for  n  >  N  and  in  the 
limit  of  large  n(x),  Pj^(x)  =  1. 


3.  The  Probability  of  "Maximum-Pack" 

It  is  clear  from  the  physical  arguments  and  from  Equation  (6)  that 


there  is  a  non- zero  probability  PR(y)  of  counting  the  maximum  number, 

N,  of  vehicles  which  can  be  packed  into  an  interval  y(NxQ<  y  <  (N+1)xq). 

The  probability  pn(NxQ)  is  identically  zero  since  counting  begins  with 

the  passing  of  a  vehicle.  However,  if  we  look  at  an  interval  of  length 

♦ 

Nxq  with  randomly  chosen  origin,  the  probability  q^NxJ  of  counting 
the  maximum  number  in  Nxq  is  the  convolution  of  the  probability  that 
the  first  one  appears  at  s  and  the  probability  Pn.  i(Nx0-s)  that  (N- 1) 
vehicles  appear  in  the  remainder  of  the  interval. 


qN1Nxo)=(xo+  ,A"1)'1  J 


pN-  i(Nxo'  s)  y 


atx)dx  ds 


(14a) 


,  ,  ,  *o  y[N  -  1;  p(x  -  s)] 

qnCNxo? =(xo+ p  )_  j  — (n  -  in -  d8  (14b) 

0 


(x)  corresponds  to  U  (x)  in  Reference  &  ;  i.  e.  ,  a  randomly  chosen 


counting  origin. 
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(14c) 


00  00 


l  V  V  e'^Xo  (px  )l 

*  <‘ +  ■“„>  J  I  - 

j=N-l  i=j+l 


l  ! 


In  two  limiting  cases  we  note  that 


Lim  qN(NxQ)  =  1 
p  - oo 


(15a) 


and 


q 


l(xo) 


|1X 

n  o 

1  +  px 
r  o 


(15b) 


The  first  of  these  two  equations  restates  the  rather  obvious  physical  fact 

that  in  the  case  of  regularly  spaced  traffic  the  probability  of  maximum  pack 

equals  unity.  The  second  expression  can  also  be  derived  by  noting  that  the 

probability  of  finding  one  vehicle  in  is  equal  to  the  fraction  of  the  average 

inter -vehicle  spacing  occupied  by  an  interval  of  length  xq  or  xq(xo  +  p  )  = 

=  ux  ( 1  +  ux  f  1  .  Figure  3  is  a  plot  of  the  probability  of  maximum  pack, 
o  o 

qN(NxQ).  as  a  function  of  f*Xo,  the  ratio  of  the  minimum  headway  to  the 
average,  1/p,  of  the  exponential  portion  of  the  inter-event  distribution  of 
Equation  (1). 
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SUMMARY 


Figure  4  is  a  plot  of  two  experimental  curves  of  the  cumulative,  AQ(t), 
of  the  probability  distribution,  a(t),  of  Equation  (1). 

Figure  4(a)  is  a  plot  of  the  fraction  of  cars  east-bound  on  Upper  Market 
Street  having  headways  greater  than  t  seconds.  These  experimental  data 
were  obtained  from  the  living  room  of  251  Upper  Terrace,  San  Francisco  17, 
in  September  I960. 

Figure  4(b)  is  a  plot  of  the  cumulative  of  the  distribution  of  inter-plane 
arrival  times  over  Runway  13L  of  the  Midway,  Chicago  Airport.  VFR  con¬ 
ditions  were  observed  on  the  6th  of  February  1959  when  these  data  were 
taken. 

With  the  kind  permission  of  P.  K.  Kinzbruner  I  am  reprinting  one  of  his 
experimental  curves  measuring  the  fraction  of  cars  having  inter -vehicle 
headways  greater  than  a  given  size.  These  data  were  obtained  in  March  and 
April  of  I960  in  Boston,  Massachusetts.  Further  details  of  the  experiments 
are  available  in  his  thesis.^  ^ 

The-  fact  that  the  expressions  for  the  counting  distributions  also  support 
other  experiments  not  reported  here  suggests  that  the  effect  of  minimum  inter 
vehicle  spacings  and  headways  may  be  important  in  medium  and  high-density 
traffic  flows. 
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NUMBER 


0  1  2  3  4  5  6  7  8  9  10  11  12  13  14 


t(seconds) 


Fig.  5  --  NUMBER  OF  INTER-ARRIVAL  TIMES  GREATER  THAN  t  SECONDS. 
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